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ABSTRACT 


197 Au between 5 MeV and 40 MeV 


Inelastic electron scattering from 
excitation energy revealed giant multipole resonances. Angular distri- 
bution studies with incident 90 MeV electrons show evidence of previ- 
ously unreported giant monopole states. Reduced matrix elements have 
been extracted and multipolarity assignments have been made. Seven 
States have been observed at excitations of 7.3 (undetermined magnetic), 


9.2 (EO), 10.8 (E2), 14.0 (El), 18.0 (undetermined), 22.5 (EZ), and 
33.5 MeV (EO or Be) 
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I. INFRODUCTION 


Giant resonance investigations of heavy nuclei have been actively 
pursued at the Naval Postgraduate School electron linear accelerator 
since 1973. Initial studies by Warshawsky and Webber [Ref. 1] demonstrated 
the feasabiltiy of measuring inelastic electron spectra in the continuum 
region of 197,y, They adapted analysis techniques developed at Darmstadt 
and measured cross sections for several known resonances. Warshawsky 
and Webber further suggested the existence of previously unreported 
resonances at about 9 MeV and 23 MeV excitation energy, but their anal- 
ysis did not conclusively identify the multipolarities. 

The work described in this thesis completes the study of 197 Ay 
discussed above and extends the excitation energy range to 40 MeV. 

Seven giant multipole resonances have been observed. Besides the pre- 
viously reported 7.3, 10.8, and 14.0 states [Refs. 1-4, 18 and 19], the 
existence of resonances at 9.2 and 22.5 MeV are required to explain the 
observed electron spectra. 

Experimental values of the inelastic scattering form factors were 
measured at scattering angles of 60°, 75°, 90°, 105°, and 120° with 
incident electrons of 90 MeV energy. Transition multipolarities have 
been deduced and reduced transition probabilities have been extracted. 
Newly identified are resonances at excitation energies of 9.2 MeV 
(isoscalar EO), 22.5 MeV (isovector E2), and 33.5 MeV (isovector EO 
or E2). Previously reported multipolarities for the 10.8 MeV (isoscalar 
£2) and 14.0 MeV (isovector F1) [Ref. 18] resonances have been confirmed. 
The analysis contained here is unable to uniquely assign transition 


multipolarities to the 7.3 MeV (magnetic) and 18.0 MeV structures. 





Pe EORY 


A. ELECTRON SCATTERING 
The Mott cross section for an extremely relativistic electron 


which scatters elastically from a nucleus of point charge Ze is given 





by [Ref. 6] 
2 2 2 
aie: ( ze | Cos z 
io = === . 4 II-] 


Nott 


The electron is assumed to have spin while the nucleus has neither 
spin nor magmetic moment. Nuclear charge structure in a nucleus of 


finite extent modifies the scattering cross section to [Ref. 6] 


i LEA 





cer 
( 2 \ = (32 | || awe Par 
dz / ‘ 
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volume 


My 


e (7) is the charge density within the nucleus as a function of T, the 
radius vector from the center of the nucleus, and @ is the momentum 
transfer vector. The Born approximation has been assumed. The temn 


describing the nuclear structure is called the square of the form factor, 


S. iS ; 
{- 2>) Moagure 
| CF (5 , ila 





For inelastic scattering, the cross section may also be written as a 
modification to the Mott scattering cross section, but the form is 
more complex, with both transverse and longitudinal contributions 
which are described elsewhere {[Ref. 7]. For inelastic longitudinal 


Searrering the form factor 15 related to the matrix element of the charge 
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operator evaluated between initial and final nuclear states. Although 
transverse components require modification of the simple expression 
II-2, it is conventional to call the experimental form factor squared 
the ratio of measured to Mott cross sections as in II-3. This value 
1s compared to theoretical calculations based on particular nuclear 
models. 

For inelastic scattering, the Born approximation may be used to 
compute the theoretical form factor based on an assumed nuclear transi- 
tion charge distribution. First attempts at a solution were made by 
describing the perturbing fields in terms of a multipole expansion, 
and assuming the incident and scattered electron wave functions to be 
plane waves. This so called plane wave Born approximation (PWBA) is 
satisfactory for light nuclei, but becomes less accurate as the number 
Sbeprocens iwie Mucieus increases. For heavy nuclei, the Coulomb 
field surrounding the nucleus distorts both the incident and scattered 
electron wave functions sufficiently that a plane wave solution no 
longer gives the correct results. Consequently, the plane wave solutions 
of the PWBA have to be replaced by phase shifted spherical waves. This 
method, termed the distorted wave Born approximation (DWBA), in practice 
requires computer computations to extract the form factors ([Ref. 8]. 
Form factors thus calculated are plotted as a function of scattering 
angle for each desired electric and magnetic transition. Comparison 
of these plots with experimental results allows for the assignment of 


transition modes to the observed resonances. 


B. GIANT RESONANCES 
The term ''giant resonance" has been used until recently to describe 


the giant electric dipcle resonance (GDR). This resonance appears as a 
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well known feature in the nuclear excitation spectra produced by such 
photonuclear processes as photodisintegration, photon absorption, and 
inverse reactions such as the (p,g& ) reaction. The presence of the GDR 
is characterized by a giant bump in photonuclear spectra at excitation 
energies from 10 to 25 MeV. It occurs throughout the periodic table 
and, according to the Goldhaber-Teller and to the Steinwedel-Jensen 
model, appears at excitation energies that are proportional to the 
inverse nuclear radius, ES. os soa 1/5 [Ref. 9]. 

This conspicuous resonance appears also in inelastic electron 
scattering spectra. The advantage of electron scattering experiments 
over the photonuclear processes for exploring the nature of the giant 
resonances lies in the variable momentum transfer, q. In the case of 
the photonuclear processes, ne momentum transfer is determined uniqucly 
by the excitation energy. The ziant resonance in the electron scattering 
cross section can be measured at different values of q depending on the 


scattering angle,@, and the incident electron energy, E:; the dependence 


1? 


being given by 


NTO) 


2 . 46.E sin* 
vie ae i + 
II-4 


where Er is the final electron energy. In addition, at values of q 
higher than is available from photonuclear work, excitation modes 
other than El become observable. Consequently, the term "giant 
resonance" now refers to all multipole excitation modes which exhibit 
Strongly coherent properties. 

Early attempts to explain the giant resonances were based mainly 
on collective models. The initial Goldhaber-Teller model assumed a 
rigid displacement of the protons as a whole against the neutrons 


as a whole, thus producing a large collective dipole moment. The 
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Generalized Goldhaber-Teller model assumes that nuclear matter consists 
of four interpenetrating fluids: neutrons with spin up (nf), protons 
with spin up (p#), neutrons with spin down (nJ}), and protons with sri 
down (p¥). Any two of these fluids may oscillate 180° out of phase 
with the other two, giving rise to four possible collective modes. For 
example, nucleons with spin up could oscillate against nucleons with 
spin down producing a magnetic resonance’ ([Ref. 5}. 

The Steinwedel-Jensen model or hydrodynamic model is based on the 
oscillation of neutron and proton fluids within a rigid nuclear surface. 
An extension of this model is the dynamic collective model ([Ref. 9]. 
The dynamic collective model replaces the rigid nuclear surface of the 
Steinwedel-Jensen model by a surface that may vibrate. These surface 
vibrations, however, are much slower than the internal giant dipole 
oscillations sa that the approximation can be made that the dipole 
oscillations take place in a deformed nucleus. 

A different model of the dipole giant resonances,used mainly for 
light, closed-shell nuclei, is based on the shell model of the nucleus. 
This particle-hole (ph) or Brown model arises from the observation by 
Halkanson (Ref. 10] that a lp nucleon, promoted into the Zs-ld shell, 
will produce dipole strengths large enough to explain the observed giant 
resonance cross sections. The shell model description [Ref. 10] gave 
the strength of the absorption correctly, but the resonance energy was 
incorrect, and the effect of the particle-hole interaction was noted 
byaBrown (Ret. 11]. 

It has been shown [Ref. 9] that such a shell model is separable 
in the relative co-ordinate of all the protons versus all the neutrons; 
and that, upon absorption of a dipole photon, a particular linear con- 


bination of single-particle excitations is produced which corresponds 
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to a motion of all the protons against all the neutrons; i.e., a 
Goldhaber-Teller collective vibration. Hence, the shell model and col- 
lective models are equivalent in their description of collective states. 

Giants resonances, then, may be described as strongly collective 
modes of excitation in which a considerable number of nucleons take 
part. In the last several years, much progress has been made in the 
investigation of these states using models such as those described 
above. Several modes of excitation other than the GDR have been identi- 
fied, including the giant electric quadrupole resonance (GQR) and the 
glant magnetic dipole resonance (GMDR). A short summary of some of 
these works follow. 

Early electron scattering investigation of the giant resonance 
structure was carried out by Isabelle and Bishop on 10) [Ref. 7]. 
By studying the angular dependence of cross sections, the giant dipole 
was observed. These experiments demonstrated the feasibility of systematic 
studies of the level schemes and transitions in light nuclei. 

emapwWeni et. al. [ReraiZ] have carried out extensive investigations 
of the fine structure in the giant resonance region of oe By making 
measurements at both backward and forward angles while keeping the 
momentum transfer constant, they were able to extract the longitudinal 
and transverse form factors as functions of excitation energy. The 
resolution of the fine structure in this manner revealed seven peaks 
associated with the longitudinal form factors and nine with the trans- 
verse. The dipole strength is distributed among a number of states to 
produce this fine structure in the giant resonance region. 

An investigation of the isotopic characteristics of the giant dip :c 


Eesomance has been carried outeby Carlos et. al. [Ref. 13), The giant 
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i 146 148 


dipole resonances for the Neodymium isotopes teu to Nd, Nd, and 
150vq were studied with a monochromatic photon beam by measuring the 
partial photonuclear cross sections O(¥,n) @(¥, pn) ando(¥,2n). 
As the number of neutrons in the nucleus increased, the dipole resonance 
began to broaden. For Nils the resonance actually split into two 
distinct peaks. An interpretation is given by the dynamic collective 
model in which there is a coupling between dipole oscillations and 
surface vibrations. The splitting of the resonance in Na characterizes 
a deformed nucleus, and the experimental cross section was closely 
reproduced by superposition of two separate Lorentz lines whose centers 
are separated by four MeV. In the framework of the hydrodynamic model, 
one can assume for an axially deformed nucleus two modes of excitation 
at differing energies, equivalent to oscillations in the direction of 
the two axes. 

Pitthan and Walcher [Ref. 14] did (e,e') investigations with Cc, La, 
and Pr. Al] three elements show distinct resonances at excitation 
energies of 9, 12, and 15 MeV. The resonance at 15 MeV was identified 
as the El giant resonance based on correlation of the excitation energy 
and strength with the GDR seen in photo-neutron experiments. The 12 
MeV resonance exhibited longitudinal properties similar to the E] 
but did not appear in (@,n) studies. It was assigned an E2 or EO 
excitation mode based on comparison to form factors developed from the 
hydrodynamic model. The 9 MeV resonance was found to be transverse 
and most likely an Ml transition. 

Later, Lewis and Bertrand (Ref. 15] observed that the (p,p') spectra 

ee) She | 120., 


, and Ba consistently produced enhancements at 


excitation energies slightly lower than those of the known giant 
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dipole resonance. Since the resonance energies of the giant dipole had 
been previously well measured, they concluded that it was unlikely that 
the giant dipole excitation was responsible for most of the observed 
enhancement. This apparent resonance shift was due to the presence of 
a giant isoscaler-quadrupole vibration located at EL = 63 a3 MeN. 
In forward angle (e,e') experiments, Fukuda and Torizuka observed 
[Ref. 16] similar El] and E2 seem in ise The usual dipole reso- 
nance was seen at 16.65 MeV and additional resonances were seen at 14 MeV 
and 28 MeV. Based on a prediction of the hydrodynamic model that a 
quadrupole giant resonance should appear at a position 1.6 times the 
peak energy of the usual giant dipole resonance [Ref. 17], the 28 MeV 
resonance was assigned an E2 transition. An assignment of either E0 


or E2 was made for the 14 MeV peak. Based on similar results for 22 Eee 


116 08 


o) 


Z ; : 
Sn, and Po (Rez. 15] these new azant resonences seem to be wiiversal 
f oo : . : -] =I 
features existing at excitation energies of 65 A / Mey and 130 A ley 


1/3 


in all nuclei. Table I contains a listing of the A> dependence of 


the resonances identified in this work and a list of references where 
others have scen or predicted corresponding resonances in 197qu and 
other nuclei. 

TABLE 1f 


= 


A Dependence of Eos Resonances 


E, @leV) EB (A0/°-MeV) References 


53 t {Refs. 21 § 22] Ce, La, Pr [Ref. 23] 
18.0 105 Pb [Refs. 24 §& 26] 


eee 3 oy) Miners Seer [Ref 24] 
25.15 105 eeciec. |e PDT Ret. 26] 


t = Theoretical Prediction 
16 





C. TRANSITION PROBABILITIES AND SUM RULES 


The transverse-electric and Uaeneere form factors are given by [Ref. 27] 


| Fe yl (4) = = ( a E= Bla! tJ, Jd.) In'S 


The continuity equation provides similar relationships for the longitudinal 
form factor. The coefficients B are the reduced nuclear transition prob- 
abilities and, as such, contain all of the nuclear physics of the inelastic 
scattering problem. The theoretical form factors are calculated by the 
GBROW computer program as described by Ziegler ([Ref. 27]. 

To determine whether or not the resonances observed are collective 
phenomena, comparison of experimental to single-particle reduced tran- 
sition probabilities may be made. For a heavy nucleus, one would expect 
a strensth greater than several times the single-particle value. Further- 
more, a giant resonance should exhaust an appreciable fraction of the 
appropriate sum rule. If an observed resonance greatly exceeds the sun 
rule under the assumption of a particular transition multipolarity, it 
is unlikely that the multipolarity assignment is correct. 

The single-particle reduced transition strengths are given by the 


Weisskopf unit [Ref. 28], 


_ & Quy _3- R | leo 
BEL). 7 Sir ie 


and 





B(mu) = e(2U+1 S =( 
P 


Ww L+3 


ey 





-1/3 197 
where R= 1.2 A = 6.98 fm for Au. Table II provides values of 


the single-particle transition strengths for ee 

Another, to some extent model-independent, evaluation of the transition 
Strength of EL modes may be achieved by expressing the strength relative 
to the energy-weighted sum rule (EWSR). Sum rules for magnetic transitions 
are more dependent on the structure of the nucleus in question and, there- 
fore, no comparisons to these have been made here. States at excitation 
energies less than 15 MeV appear to be AT = 0 except for the EL mode 
which is AT = 1. The sum rule for isoscalar (AT = 0) excitation modes 
with L?1l is given by [Ref. 29] 

2,2 2b-2 
S (EL AT=0)= S16 £8 (ELF) = Fe Labey WER PD II-8 


Tr AM 


il 


The isovector (AT = 1) sum rule for L>1 is related to the isoscalar 


sum rule by 


N 
S (EL aT=1) S (EL 4T=°) > . II-9 


The corresponding sum rule for an isoscalar monopole (E0) excitation is 


fRef. 30] 


S (EO) = oe | Mel 


ff 


hz 2 p? 
hE CR 11-10 
eae 


where Mes is the monopole matrix element. The strength of El resonances 


can be expressed relative to the sum rule given by Warburton [Ref. 31]: 


ele (3 Ne 
S(el) = Zmp lar] A aah 


Es) , - 
Table III gives values of these sums for Au using weighted radii of 


the ground state charre distribution. 
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TABLE II 


7 

Weisskopf Units for MY Au 

Gell, B(ML) * 
L et p2 
A 6.6 .06 
2 3, 35x10" 3.1 

4 2 

3 1.63x10 1.44x10 





* Units are fm 


TABLE IIT 


Sum of EWSR For 197 Ay 


EL AT Sum Units 
9.2x104 MeV-fin' 
7.04x10° e*-McV-£m2 
7.4x10° e*-Mev- fn! 
1.1x10° e*-MeV-fm4 
6.15x10° e?-MeVv-£mn® 
9.16x10° o*-MeV-fm® 
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D. QUASI-ELASTIC SCATTERING 

Electrons incident upon a nucleus may scatter elastically with 
individual nucleons within the nucleus. This process is called quasi- 
elastic scattering. If sufficient momentum is imparted to the nucleon, 
it will be knocked out of the nucleus with some kinetic energy. The 
Simultaneous energy loss by the incident electron may result in a 
measurable quasi-elastic peak.- It 1s of interest to see at which value 
Ope E this contribution appears on the spectrum of arae. EtOmsee alt 
it contributes to, and therefore must be subtracted from, the scattered 
electron spectrum in the region of the giant resonance. 

The inelastic electron-nucleus cross section has been calculated 


in three energy-loss regions with the Fermi gas model for the nucleus 


[Ref. 32]: the quasi-elastic region; the threshold pion region with its 


3-3 resonance of the proton. Pion production begins at 135 MeV and the 
3-3 resonance is 298 MeV above the elastic peak, so these regions are 
not of interest in this work. 

In the Fermi gas model, the electron scatters elastically from a 
Single nucleon in the free Fermi sea, with the recoiling nucleon required 
to scatter out of the Fermi sphere. Because of the Pauli principle, the 
nucleon cannot scatter into an already occupied state. The final nucleon 


energy 1S given by 


= 


- | ko + we}? ; II-12 
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where k is the nucleon momentum and M is the nucleon mass. The initial 


(bound) nucleon energy is 


2 2 t-* 


where the effective nucleon mass, M*, is given by [Ref. 33] 


* 


N\ = MM II-14 
1.4 

and U(k’) is the effective single-particle potential in nuclear matter. 
Thasepotential effectively shifts the electron energy loss, , to take 
into account the nuclear binding. Moni treats U(k*) as a constant, E , 
equivalent to the average nucleon interaction cnergy. The quasi-elastic 
peak energy loss will increase with q and approach the effective free 
nucleon kinetic energy, de i (2M"), at large values of q. 

The magnitude of the above shift for the cold nucleus is of interest. 
The average nucleon interaction energy, & , and the nuclear Fermi 
momentum, k.. , have been determined for nine elements from Li to 
208pp, {Ref. 34) by a least squares fit of theory to quasi-elastic peak 
data. In the cross-section formula, an energy-conserving delta function 
appears, given by 


8 (Ww + _ -€- oS) BIG 


) 
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which determines LJ . The nucleon momentum is approximated by the Fermi 
momentum and is roughly constant at 260-265 MeV/c for heavy nuclei from 
nickel through lead. Cc fee 2 oll (Ref. 34]. The energy-conserving 


delta function requires that 


car ya 
ee ee II-16 
NN 2M 
For the smallest momentum transfer of these experiments, W = 73 MeV. 


The threshold value is W= Crus 43 MeV. Thus explicit corrections 
for the quasi-elastic contribution need not be made for the spectra of 


this work. 
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III. DATA ANALYSIS 


A. DATA REDUCTION 

Five experiments were performed with an incident beam energy of 
90 MeV and at scattering angles of 60°, 75°, 90°, 105°, and 120°. 
The experimental information is given in Table IV. In general the data 
analysis procedure described by Warshawsky and Webber [Ref. 1] was 
followed throughout this work. The ees choice of resonance energies 
was determined from known resonances [Refs. 1, 2, 3, and 4] and by visuai 
inspection. Previously determined resonances occur at ] See Oe 
and 14.0 MeV. Visual inspection of the background subtracted data 
revealed wide structures at 18.0, 22.5, and 33.5 MeV. 

ihe fitted anelastic spectra extendedsirom 4 to 40,.Mev excepe 
the 120° spectrum which extends oniy to 51.5 MeV. In some cases the 
data extends to 50 MeV but no structures were observed at this high 
excitation energy. Since the 120° data has different boundary conditions 
for matching the background at high excitation energies, only cross 
sections for resonances below 20 MeV were extracted at this angle. 

Breit-Wigner resonance forms were assumed for each resonance and 
computer code NAW (see appendix) simultaneously fits the entire inelastic 
spectrum including the radiation tail and background. The new resonance 
positions and widths were manually varied while the NAW code varied the 
peak strengths and background parameters until a best fit was obtaincd. 

The criteria used to detemnine a reasonable fit were as follows: 

1. The data and calculated spectrum should coincide visually. 


Ze x? per degree of freedom should be less than one. The data is 


not strictly stutistica] because the detector momentun interval is larger 


Zs 








than the momentum increment of the spectrometer field and hence correla- 
tions exist between energy bins. 

3. The difference between the coca Xx : for successive computer 
iterations wasat’ 2 0.5. Typically there were 300 degrees of freedom 
and a change of 0.5 in the total chi squared resulted in 17 x one per 
degree of freedom. 

4. All observed resonances should consistently fit each of the 
five spectra. Table V lists the positions and widths of the resonances 
required. 

Figures 1 - 5 present the experimental data corrected for spectrometer 
dispersion effects with the fitted total background and the individual 
resonances drawn. Figures 6 - 10 present these spectra with the total 
background subtracted, revealing the resonance structures of Table V. 

The areasunder the inelastic resonances were calculated from the 
Breit-Wigner resonance shapes and the ratio of each inelastic resonance 
area to the area under the elastic peak was calculated (A;/A,). The 
inelastic form factors squared are obtained by multiplying A. 1Ke by 
the square of the elastic scattering form factors (oP Table VII 
presents the experimental ratios of the inelastic to elastic areas and 
the square of the inelastic form factors. The elastic form factors 
were Calculated using the program of Rawitscher and Fischer [Ref. 35]. 


The charge distribution parameters were c = 6.38 fm and t = 2.32 fm. 
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TABLE IV 


Experimental Conditions 


Target Incident 
Thickness Energy 


048 g/cm 90.32 MeV 


096 g/cm 90.24 MeV 


096 g/cm” 90.21 MeV 
144 g/cm? 89.98 MeV 
240 geome $9.60 MeV 





A TXT 


TREE v 


Observed Resonances 





ZS 
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FIGURE 1 
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TABLE VI 


Scale Factors for Converting 
Arbitrary Units to Expcrimental Units 


Angle Scale Factor MeV-sgr / Arb. units 


9.9 x1074 


1.89x1074 


2.99x107> 


1.53x107> 


1.47x10"> 
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TABLE VII 


Area Ratios and Inelastic Form Factors Squared 


57 


Total 


% Error 





— 


-30x1074 


-4 
80x10 


07x10" 


-40x107> 


ecilecoed 








TABLE VII Continued 





7.60x10 > 







1.12x107¢ 
5 












4,.42x10 





$.46x10 > 


3.18x10°- 





1.82x10°> 


5 











2.70x10— 
3 





ier caKas 


6.91x10°° 


2 









3.41x10- 


4.21x10"> 


~3 
7.83x10 


7.17x10° > 


5 








8.60x10_ 
















-3 =A 
7 O60 9.09x10 
9.80x10 > 4.35x10°" 
9.90x10°> 2.48x10°" 

oy) -4 
1.5xl0 3.2x10 
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B, ERROR ANALYSIS 

The procedure outlined in the previous section provides values of 
the ratios of A,/A, but the uncertainties attached to these ratios 
requires additional considerations. The program NAW computes only 
Statistical errors. There are other errors, both systematic and random 
which are not strictly statistical. Some of these are the errors result- 
ing from the instrumental effects and are outlined by Bernard and 
Traverso [Ref. 36]. 

To obtain a more realistic estimate of the errors in the ratios 
A./A. the non-statistical uncertainties are assumed to be those which 
arise from the uncertainties in locating the resonance positions and 
determining their respective widths. These arise from two sources: 
the overall instrumental uncertainty which 1s estimated to be 0.15 MeV, 
and the frecdom allewed in the resonance peak positions and widths in 
fitting these spectra. The total uncertainty of the energy in this work 
is the sum of the fitting possibilities and instrumental errors. 

The Breit-Wigner shape is then analysed to see how these variations 


in energy effect the area under the resonances. The Breit-Wigner shape 


1s given by 


AOCELE TM) A 
ei |S (E-E,+(%5Y ? III-1 


Gh), 





where EA is the resonance energy, \ is the resonance width, A is the 
Strength parameter, and E is the excitation energy. The resonance cross 
section is given by the integral of the Breit-Wigner form 


-|, Ep 
a / 
O~(E, \') = — + 2° = (7) tie? 





The variations in position and width introduce an uncertainty in 


the resonance cross section, 





oo So 
o- ~ |——jAE, +] AV : 
Bye, .t ) Se R | \P | III-3 
The absolute values of the partial derivatives are assumed since the 
peak position and width uncertainties may be correlated. 


The fractional error in the resonance _ section 1S 


5 Wy 2 Tan Ct) me, Le ¢ 
eae [rie |aee + | Ae 


Ae 2. Tart (Ee/p) 








Lies 


This uncertainty is then added to the statistical error to provide 


the total error. Table VIII lists the uncertainties from the two sources. 
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TABLE VIII - 


Fractional Error Contributions 


a a 





o_ 


Limit 27 
Upper oly 
Linit 

TLD sill, 
25 OS 
» 34 

Upper 

Limit 
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IV. DISCUSSION 


Pe COLLECTIVITY 

Data analysis techniques explained in Chapter III were used to sub- 
tract the background and extract the elastic and inelastic spectra 
of the scattered electrons. A least squares fit was used to evaluate 
these spectra and resulted in the identification of seven giant 
Messrances at excitation energies of 7.35, 9.Z, 10.8, 14.0, 18.0, 
22.5, and 33.5 MeV. These resonances have strengths corresponding to an 
appreciable fraction of the appropriate sum rule as exhibited in Table IX. 
Additional support to reas collective nature is given in Table X where 
it 1s seen that the transition strengths are greater than several Weis- 


skopf units. True giant resonances are therefore being observed. 


B. MULTIPOLARITY ASSIGNMENTS 
7.3 MeV Resonance 

The multipolarity assignment of the resonance at 7.3 MeV is un- 
eemedine Buskirk et. al. [Ref. 18] and Lone et. al. (Ref. 19] have 
observed this structure in backward ancle (e,e') experiments and indicate 
a transverse excitation. A comparison of the data with calculations using 
GBROW [Ref. 27] and a code by Drechsel |[Ref. 37] is incapable of deter- 
mining a definate multipolarity assignment. Among possible reasons for 
this situation are a multipolarity admixture, existence of unknown 
longitudinal resonances influencing the forward anple data, and experi- 
mental difficulties encountered by the high counting rate at forward 


scattering angles. 


A2 





9.2 MeV Resonance 


iy 


The analysis of the 9.2 MeV resonance (53 A. 2 favors an EO 
multipolarity assignment. This resonance exhausts 38% and 18% respect- 
ively of the EO and E2 energy-weighted stm rules and its strength 
corresponds to five Weisskopf units when an E2 assignment is made. 
However, the angular dependence of the experimental form factors favors 
EO (Fig. 12). The result would be in agreement with the isoscalar 
giant monopole resonance predicted at 56 Ay [Ref. 21] and at 
58 al/s MeV [Ref. 22]. Resonances exhibiting the characteristics of 
EQ or E2 transitions were also seen at Ex = 53 ie Sine te N= ae 
maecley, (Ref. 25] and in ac thy [Ref. 24]. An El multipolarity assignment 
is ruled out by comparison with DWBA form factor calculations and further- 
more by the fact that no resonance has been seen in © - absorption 
measurements [Ref. @] at this excitation energy. 
10.8 MeV Resonance 

The resonance at 10.8 MeV (63 a l/3) was seen previously in beg 
Refs. 3 and 18] and in many other nuclei (sec e.a. [Ref. 5]). Generally 
an assignment of E2 was favored in these cases. In this experunent a 
transition strensth of 14 Weisskopf units is obtained assuming an 
electric quadrupole excitation. Comparison of the experimental form 
factors to the angular distribution calculated from the DWBA form factors 
is. inconclusive in distinguishing between an EO and E2 modes (Fig. 13). 
However, the experimental strength corresponds to 113% and 68% of the 
EO and £2 sum rules, respectively. Together with our results on the 
9.2 MeV resonance, the E2 assignment is therefore favored. This resonance 


is thought to be the isoscalar quadrupole state predicted by the 


collective vibration model of Bohr and Mottelson [Ref. 38]. 





14.0 MeV Resonance 

The 14 MeV (81 A71/3) giant dipole resonance was measured with 
monochromatic &- experiments [Ref. 25] and has also been seen in other 
{e,e') experiments [Refs. 3 and 18]. For this experiment, form factor 
comparisons to DWBA calculations show reasonable agreement at forward 
angles, but exhibit deviations at backward angles. An attempt was made 
to further explore the nature of the giant dipole transition by assuming 
both volume and surface oscillations separately in the DWBA calculations 
(Fig. 14) [Ref. 23]. Transition strengths of six and eight Weisskopf 
units and 84% and 104% of the sum rule are obtained for the volume and 
surface oscillation modes respectively. However, neither model adequately 
describes the observed data. As suggested by Pitthan [Ref. 39], the 
discrepancy at backward anglesmight be explained by an electric spin flip 
contribution to the El resonance [Ref. 40]. Using only the forward angle 
data to extract the reduced matrix elements, the volume oscillation 
strength exhausts the sum rule (98%), whereas the surface oscillation 
exceeds the sum rule. Thus a volume oscillation is favored. 

18.0 MeV Resonance 

A resonance of undetermined multipolarity is observed at 18.0 MeV 

(105 alls, excitation energy. Nagao and Torizuka [Ref. 24] have observed 


BUS but the state in 127au 


E3 strength at the corresponding energy in 
does not conform to this multipolarity. E(O-3) and M (1-2) calculations 
‘have been compared to the data (Fig. 15), but no reasonable agreement 


is found. It appears that the cross-section at this excitation energy 


Some nseamixture of vaaious transitions. 
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22.5 MeV Resonance 

Comparison of the experimental results for the 22.5 MeV (1307 1/3) 
resonance with DWBA calculations (Fig. 16) gives support to an E2 
assignment. At forward angles the data is indistinquishable between 
the EO and E2 transitions, but the 105° measurement is incompatable 
with the EO assignment. The Sendai group has reported [Ref. 3] an 
electric monopole or quadrupole at this excitation energy in 17a and 
also at the corresponding energy in 208pp [Ref. 24]. Other multipolarity 
assignments are unsuccessful in predicting an angular distribution in 
eeyeenent with the data. Ihe 22:5 MeV state 1s thought to be the 
isovector quadrupole state predicted by the collective vibration model 
of Bohr and Mottelson [Ref. 38]. 

33.5 MeV Resonance 

A new resonance, previously unreported in vere is observed at 
33.5 MeV (194 A“l/5), This corresponds roughly to the isovector giant 
monopole resonance predicted by Suzuki [Ref. 22] through calculations 
based on the collective vibration model of Bohr and Mottelson. Both 
the EO and E2 assignments exceed their respective energy-wcighted 
sum rule prediction. However, our line shape fit, while very sensitive 
to any structure in the spectrum, may tend to over-estimate cross- 
sections for very wide resonances. For this reason, the B-values may be 
too large. Comparison of the experimental form factors to the DWBA 
calculations allow both EO and E2 multipolarities to be possible. A 


208 
Pb at a corresponding energy of 


Similar new resonance was observed in 
32.5 MeV [Ref. 26] and tentatively given an EO or E2 assignment. It is 
toebeeemonasized that consistent fits for all the spectra (97 au) 

over the entire range from 5 MeV to 40 MeV excitation energy were only 
possible with the inclusion of this resonance at 33.5 MeV. 
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TABLE X 


Transitions Strengths in Weisskopf Units 
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TABLE X Continued 
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V. CONCLUSIONS 


Giant multipole resonances in 1 were studied with inelastic 
electron scattering. Resonances at 7.3, 9.2, 10.8, 14.0, 18.0, 22.5, 
and 33.5 MeV excitation energy were observed. Of these, the 9.2, 18.0, 
22.5, and 33.5 MeV states have not previously been reported in the open 
literature. Multipolarity.assignments of these states were made as 
was discussed in section IV-B. 

The 14.0 MeV state is identified as the much studied giant dipole 
resonance, see e.g. [Refs. 3, 5, § 25]. Our analysis suggests that this 
El state may be described by a volume oscillation of the charge distribution. 
Furthermore, the data 1s compatible with a suspected transverse spin flip 
contribution to the El cross section. 

The angular distributions of the cross-sections for the 10.8 MeV 
and 22.5 MeV resonances are compatible with previous E2 multipolarity 
assignments [Ref. 3]. In addition, this investigation adds evidence 
for the exclusion of a monopole assignment for the 10.8 MeV state. 

_ 4197 208 

A state was found at 18 MeV in FAV Pa Bo Pb a resonance was 
seen [Ref. 26] and cross-section analysis of another (e,e') experiment 
found E3 strength between 16 MeV and 22 MeV [Ref. 24]. The analysis 
of the experimental results for this energy region in re is not 
capable of assigning a definite transition multipolarity. The data 
suggests admixtures of longitudinal and transverse contributions to the 
GeGeS-SCCLION . 

The 7.3 MeV transition is known to be magnetic. However, this 
investigation, like others [Ref. 19], was not able to determine the 


transition mode. 
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Theoretical Meera con concerning the existence and location of 
the giant monopole resonance is found in the literature [Refs. 21 and 22]. 
New resonances, discovered in this experiment at 9.2 MeV and 33.5 MeV, 
have angular distributions in agreement with monopole assignments. Their 
excitation energies are comparable to theoretical predictions [Ref. 22] 
of the location of the isoscalar and isovector monopole states. For 
the lower state, an HO transition is favored over an [2 assignment. 
However, for the 33.5 MeV state, the transition multipolarity is 
inconclusive. The excitation energies and transition multipolarities 


are presented in Table XL. 
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TABLE XI 


Multipolarity Assignments 


Excitation Transition 
Energy (MeV) Multipolarity Comments 


undetermined probable magnetic, 
mixture 


EO Isoscalar 
E2 Isoscalar 


el Isovector, volume 
oscillation 


undetermined 


E2 Isovector 


ae 


Isovector 
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APPENDIX A 
PROGRAM NAW, SPECTRUM FITTING CODE 


A. GENERAL FEATURES 

Naw is a Fortran program developed at Darmstadt [Refs. 41 § 43] 
which uses analytical expressions to approximate the inelastic and 
elastic line shapes and radiation tail to best fit the experimental 
spectra. The elastic line shape is assumed to be two half gaussians of 
differing half widths matched to a hyperbolic shape approximating the 
radiation tail. Work done by the experimenters at the Institut Fur 
Technische Kernphysik at Darmstadt have found this parameterization to 
be an adequate estimate of the shape of the elastic peak (see Figure 18) 


ese ee 
gaussian 


gaussian ale 


hyperbolic 





Hae 1: Low E 


relative zero 
peak position 


ASSUMED ELASTIC SPECTRUM SHAPE 
FICURE SES 


The elastic radiation tail results from processes occurring during 
Blastic scattering that preduge lewer enerpyselectrons. The radiation 


tail is calculated from the theory [{Ref. 45] at several designated cnergies 
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and these theoritical values are fitted with a polynomial of the fourth 
degree. The total background, which includes the radiation tail, is then 
Simultaneously fit with a series of BreiaeneE line shapes representing 
the nuclear excitation structure to reproduce the experimental spectrum. 
The positions and widths of the resonances are free parameters to be 
adjusted by best fit or to be assigned from known values. The areas 
under these resonances are then calculated from which the ratios A;/Ag 
(inelastic peak area divided by the elastic peak area) are determined. 
The values A;/Ag are then multiplied by the calculated scattering cross 
section resulting in the experimental inelastic cross section for each 


of the resonances. 


B. CONTROL CARDS 

NAW uses the output from program KIK as the input data. Since 
NAW's first uSe at Te NPS IBM 360/67 [Ref.1] it has undergone various 
modifications. The present form is reproduced on the following pages, 
along with the control information. NAW accepts 200 clastic peak data 


points, 500 inelastic data points and 25 inelastic resonances per 


spectrum, 
NAW CONTROL 
ee TEXT Alphanumeric identifying remarks (20A4) 
2. NORM, PN, PNz, PN, PN, PNe PNo CELLOS) 
NORM Spectrometer energy normalization factor to change the 


Origional energy units (OU) to MeV: [MeV] = NORM* [OCU] 


ee 1«2-7, Starting values for the various parameters being 
used to fit the elastic spectrum with an analytical 
expression. 

PN, HMalf-width of a gaussian to fit the high energy side 


of the elastic peak. 





PN 


EE, THETA, D, 


Ee Energy of the elastic peak in MeV. FE is used as a 
starting value in the clastic fit. 

THETA Electron scattering angle in degrees. 

D Target thickness in g/cm. 

UNT Constant background at the high energy end of the elastic 
resonance spectrum in counts/a c. 

EEL > ELL, , GSSW) GSSW, , DSSW, E of (6F10.5) 

EEL, High energy limit for fitting the elastic peak in OEU. 

EEL, Low energy limit for fitting the clastic peak in OLU. 

GSSN Upper limit tor ealenlacin@ tiesradiation tail in CEU: 

GSSI, Lower limit for calculating the radiation tail in OEU. 

DSSW Energy step to be used in calculating the radiation 
Caren OL: 

“ref See card 7, Ly. If Eef = 0, LS ae is calculated internally. 

ares ( ZiiOeS ) 

En Higneelcctron energy lamicetor fitting thesmelastic 
spectrum 1n excitation energy (MeV). 

EF} Lowel @etron cncerey slime tor fitting the inelastic 


Half-width of a gaussian to fit the low energy side of 
the elastic peak. 


Matching point for the low energy gaussian and the 


hyperbolic segment. PN, 1s specified in McV from the 
position of the zero point and is a negative number. 


Dummy . 


PN = 0, Height of a second clastic resonance peak if 
necessary. 


Position of a second elastic resonance peak if necessary, 
only with PNg# 0. 


UNT, (4R10.5) — 


SPCCUILM IMMCXCi tatronmsencrsy in (MeV). 





i a | 2 2 oS | 
GR GR, ; GR » GRY, GR GRY (GFIN: 5) 
Regions where it iS expected that only background 
contributes to the spectrum; i.e., between resonances 
in the Pienaar spectrum. There must be at least one 
Sci OteGin. Gly input, with a maximum of three. It is 
pee ommended: to set GR uxCSsWe GRexGSsWu for giant 
resonances, energy in MeV. 


Ee oleae 22 (2203) 


L: 1 A continuous curve consisting of two gaussians and a 
segment of a hyperbolic is used for both elastic and 
inelastic fits. 

2 The clastic resonance will be fit as if L,=1 and the 
inelastic resonance will be fit using the Breit- 
Wigner formula. Each curve will be multiplied by a 
Straight line with an equation Y=1 +P. (EL- Ex (1 ae Be 
will be ae ea on card #10, Ex is ie excitation 
enerpy , and i: is the nosition of the maximum 
as specified on card #8, 

Swlnceotastile resonance will bewpit as if L)—-Z%execpt 
that NAW will determine the values P. and card #10 
will be left blank. 


L4 Parameters wnich choose the analytical expression to 
fit the background due to the radiative and experimental 
peceesses, Lois a two digit number, the units give the 
order of a polynomial e¢.¢. unit=3, P,+P.AtP. AZ, Ge eS 
are free fitted paramters. The neces chodse | the 
function to be multiplied and the areument to be used. 


< 2 : 
1£ Ly< 10 Py+PA+P A+, .4P, 9, * FSSH 
10<L,< 20 P,+PoA+P,A%+...+ FSSW 


2 2 E- sr 
20s Lo< 50 {P, aia teers +P og = exp pee | 
tao feng. THI} onan (FEY 


FSSW is the calculated radiative tail, E is the energy 

of the outgoing electron, Ee is the energy of the 

elastic peak. A-(E-E <a foro. 50,08 is read in 
Clommers pd: ine sy < 90 the analytichd- form correspond- 
ing to L5=L4- 90 1S chéosen, but A is replaced by 1/A. 


IE The number of resonances with fixed energies. If La< 0 
then the resonance will be kept fixed relative to EG), 
but 11) will be “fitted. Séescard #8, 


L4 The number of resonances with energies to be fitted. 


Note that |L.J + L, must not exceed 25. 
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10 


The number of lines widths to be fitted. 70 E.g. =3, 
means that three line widths are free and in fact are 
E (1), E.(2), and E (3). If there are more than three 
résonancts, E,(3) tirough E,(N) will all be fitted 
with the same line widths. 


For Ly e= 2 @thew first Le peaks will be fitted with 
Breit-Wigner forms, the remainders with gaussian 
shapes. For L.=0, Le= [L,(+L, 


Target geometry: 


=Q Transmission -normal thickness entered on card #3 
=1 Reflection -nomnal thickness entered on card #3 
=2Z Effective thickness entered on card #3 


Output control: 


=<0 Maximum output 
=1 Minimal output-no tables and no intermediate calculations 
z-| Limited output-no tables 


Plothne comenol : 


=) Results are not plotted. 
“7eeresuilts are plotted. 


Repetitions of NAW procedure. 


=0 No repetitions are performed. 

>0 In this case it is possible to evaluate the data with 
new parameters. ‘The new parameters will be centered 
on additional cards placed after the control cards. 


Inelastic spectrum fit iteration control. Only 
override to the default options need be centered. 


NUMBER Or Teerabions were allgenercies are fixed, 
default L442 ° 


Maximum number of iterations; default is L790. 


Minimum number of iterations; default is Lj 3=5- 
Number of iterations where width are neld constant 171 
case XY actual >3* theory. 


“(eine sracditicnetom 1aiiebe calculated usang the 
Ginsber's method [Ref. 43]. 

<0 The radiation tail will be calculated using Schiff's 
method [Ref. 44]. 





16 18 


EQ), 


VQ), 


Radiation tail fit control. ‘Only override options need 
be éntered. 


Analytical representation of the theoretically 
calculated radiation tail; default is Ly gz. 


7 Oe Z G 4 
i. SSHEP, + PoA + PA + PAA + PLA 


=2 SSW=P, + PA + PsA* + PAE + PLES 
=3 SSWP, + PoA + PE + P{EY + PSEP 
=4 SSW=P) + PE + Pah + PES + PLES 
-5 SSN=P, + PU + PU? + Psu? + Pout 


P; will be fitted internally by NAW 
A= 1/({E,] E) 

U=E - E 

Ey GSSW, + (GSSW, - GSSW, )/L1g 


<0 the radiation tail parameters P. and By will be 
specified on card #15. 


The number of parameters for the RadsaAtion tall £1¢ 
ey La <= 5); defauit is Ly 9 = 
16? default is Lig = 2, 


20 Error ellipsoid control values will not be calculated. 


See L 


= 0 


S-eetie calculation ct tne radiation tail will be 
modified by the factors entered on card #16. 


=3 Inelastic spectrum results. Card images will be 
written on file 7 and may be stored on the data cell fer 
further handling. 


IE x (2) » Se a 2) GVO): 
Resonance excitation energies in MeV; i hei those 
with variable positions first. i= jL.| a Lis LS; 

V(Z) .. 5.3) @8F10.5) 
Resonance line widths, i = [L " + Ly < 25 
Drama ta)ye= 0, thé width of the clastic line will he 


assumed. 
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eee yee 2), <5, P(i) | (SF10.5) 


He. 


AZ 


Vs. 


14. 


1s 


The P. factors when L = 2 on card #7, i = | | se) Oe eras 
1 1 5 4 


Ly heme aD, EleRAs Cal (AR LO poe 20K 321s) 

Z Atomic number. 

A Mass number. 

XRAD Radiation length in g/cm’. 

FIRA Elastic form factor for Ep, and THETA of the experiment. 

Coy. The fermi charge distribution palameters.  Since™ these 
values are used as the parameters for the trapezoidal 
charge distribution potential used in the radiation tail 
calculations in Born Approximation, these values must 
beyentencd. 

XMAS, H CASO), 

XMAS Nuclear mas cquivalent in MeV. 

hi Isotope percentage (100% = 1.00). 

A; The A. are the coefficients for the calculations of the 
ejastic form factor using a power series expansion on 
Ge 

2 : One 
Fo) = = Ayq?2 G-1) 
When A,#0 the elastic form factor will be calculated 
using the Born Approximation. 

Ur (8F10.5) 
Ujmay be used to get better starting values for the fit 
by taking the background parameters from a previous fit. 
If the U;=0, the starting values are chosen internally. 
Seome dda 7. 

Malte 2 to a. PS, E (6F10.5) 
This card is read only if igs 0. 

PJ -P5 Radiation tail parameters, see Lig for tlhe expression. 

an This value is read only when Lia 6 = 5. 
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16. XSSW, XM, XB, 


XSSW 


The data are assigned to file #4. 


This card is read only if Lo, > 0. 


(SROs S) 


Factor Multiplying the area under the elastic peak, A 
e 


Factor multiplying the M@iler term. 


factor multiplying the B 
The equation for the radiation tail is now given by 


SSW = 


MAXI 


ie HOEE 
MAXI are form 


+ XM(igller term) + XB(Brems. term)* (XSSI)A, 
factors as described by Maximon and Isabelle. 


The elastic and the inelastic 


data must be followed by one card image having the data energy step 


width entered beginning in colum 72. Repeated evaluations of the same 


set of data (i.e. 


Lig = 0) must have the following set of cards placed 


after the 14 to 16 control cards as described above; one set for each 


repetition to be performed. 


1) 
Z) 
5) 
4) 
>) 


Ly ver Lig see 
Ey Ey see 
PO) Mo EG) see 
Oy ee (9) see 
PCL) ees 5 ee) see 


cana 7 
card #5 
card #8 
card #9 


Care 0 
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